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Abstract The bar formation is still an open problem 
in modern astrophysics. In this paper we present nu- 
merical simulation performed with the aim of analyzing 
the growth of the bar instability inside stellar-gaseous 
disks, where the star formation is triggered, and a cen- 
tral black hole is present. The aim of this paper is to 
point out the impact of such a central massive black 
hole on the growth of the bar. We use N-body-SPH 
simulations of the same isolated disk-to-halo mass sys- 
tems harboring black holes with different initial masses 
and different energy feedback on the surrounding gas. 
We compare the results of these simulations with the 
one of the same disk without black hole in its center. 
We make the same comparison (disk with and without 
black hole) for a stellar disk in a fully cosmological sce- 
nario. A stellar bar, lasting 10 Gyrs, is present in all 
our simulations. 

The central black hole mass has in general a mild effect 
on the ellipticity of the bar but it is never able to de- 
stroy it. The black holes grow in different way according 
their initial mass and their feedback efficiency, the final 
values of the velocity dispersions and of the black hole 
masses are near to the phenomenological constraints. 

Keywords galaxies: spirals, structure, evolution, ha- 
los, black holes 



1 Introduction 



result by analyzing the growth of a bar instability in 
exponential disks, embedded in a Dark Matter (DM) 
halo which evolves in a cosmological context. We used 
purely stellar disks as well as disks containing both gas 
and stars. In the latter case, we studied the effect of 
only allowing the gas to radiatively cool, and we then 
examined the impact of the star formation. Our em- 
bedding technique allowed us to vary the baryon to DM 
mass ratio onto the same DM halo. In almost all the 
cases we studied, including the classically stable ones, 
we noticed the formation of a long-lasting bar. The 
onset of the bar instability, in the latter cases, appears 
to be driven by the triaxiality and by the dynamical 
evolution of the DM halo. For a few cases, cooling gas 
proved to be able to stabilize the gas, but such effect dis- 
appeared when star formation was turned on, because 
gas converts in stars before a stabilizing, dense, central 
knot can form. Therefore the problem of inhibiting the 
bar formation is still an open astrophysical problem, 
since one third at least of disk galaxies are not barred. 

The observational data indicate that massive central 
black holes exist in disk galaxies as well as ellipticals. 
Moreover black hole s masses are correlated with host 
galaxies properties; iMagorrian et al.l (|l998r ) has con- 



In a se r ies of recent papers Curir et al. ( 2006 ). Curir et aL^ 
( 2007t ). ICurir et all |2008l ) we have shown that the bar 
formation inside disk galaxies is triggered also by the 
cosmological scenario and not only by the classical in- 
stability of a self gravitating disk. We obtained such 
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eluded that the median black holes mass is 0.006 of 



the bulge mass and lKormendv fc Gebhardtl (|200lh cor- 
related the black hole mass to the bulge velocity disper- 
sion. Such a large mass concentration inside the galaxy 
could affect its structure. It is known indeed that a 
central m ass concentrations (CMC) can destroy a b ar 



Bournaud et al 



Curir et al.l (|2007l) 



(120051) ■ 

The effect of the presence of a massive black hole (BH) 
at the center of a disk galaxy could therefore mimic 
the stabilizing action of a knot of dense, cold gas which 
forms if gas is allowed to radiatively cool but not to 
form stars. 

The influence of central BHs on the dynamical evolu- 
tion of bars in disk galaxies has been examined mainly 
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bv IShen fc SellwoodI ([2004) and iHozumi fc Hernguistl 
( 20051 ) in a pure non dissipative scenario and using as 
model galaxies isolated systems. They results give very 
high values (10^'^) the minimum mass necessary for bar 
dissolution, higher than the inferred one for local spi- 
rals. 

In this paper we want to treat the same problem us- 
ing simulations of stellar- gaseous disk with star for- 
mation. We will evolve the disk model as isolated 
system and in a cosmological scenario. We will dis- 
cuss the evolution of the bar ellipticity and semi-major 
axis of the same disk endowed with different initial 
BH mass. The BH is accreting during the evolu- 
tion. Its a ccretion is regulated by a suitable feedback 
Darameter( Di Matteo et al.i ( 20031) ), which states a cou- 
pling between a percentage of the accretion radiating 
energy and the gas heating. We will also show the role 
of this feedback parameter and its interplay with the 
initial BH mass, the formation of the CMC and with 
the star formation. 

Finally, we will investigate the relation between the 
decrease of the bar ellipticity,the accretion time and the 
feedback. Moreover we will compare our simulations 
with some phenomenological issues as the Kormendy 
black hole mass-bulge velocity dispersion relation. 

The plan of the paper is the following. In Section 
2, we summarize our recipe for the initial disk + halo 
system and present the star formation recipe. In Sec- 
tion 3, we present our simulations, and in Section 4 we 
point out our results. The parameters related to the 
bars formed in the new and the old stars and to the 
global, old-|-new, stellar populations are given in this 
section. In section 5 there are phenomenological impli- 
cations and section 6 is devoted to our discussion and 
conclusions. 



2 Method 

In the cosmological simulations, we embed a gaseous 
and stellar disk inside a cosmological halo selected in a 
suitable way, as to be able to host a disk galaxy, and 
follow its evolution inside a cosmological framework: 
a ACDM model with ==0.3, nA =0.7, erg =0.9, 
h =0.7, where flrn is the total matter of the Universe, 
JIa the cosmological constant, erg the normalization of 
the power spectrum, and h the value of the Hubble 
constant in units of lOOh"^ kms~^Mpc~^. The halo is 
extracted from the cosmological simulation at a given 
redshift (we chose z = 2). The stellar -l-gaseous disk 
is embedded in the halo, in equilibrium with its grav- 
itational potential as in the isolated case; the system 
is then re-inserted in the full cosmological simulation 



and evolved. A detailed description of our method of 
producing the cosm ological scenario w here the disk is 
evolved was given in ICurir et all (l2006l) . 

As far as the isolated simulations is concerned, 
we embed a gaseous and stellar disk inside a DM 
halo, having a radial density profile of the NFW form 



Navarro et al. (|l996l) and follow its evolution for 10 
The halo has a mass of 1.36 x IO^^H^^Mq and 
4.5. The disk is set in equilib- 



Gyrs. 

initial concentration c= 
rium with the gravitational potential of the DM halo, 
but not vice-versa; as a consequence, there is a transient 
period in which the halo responds to the disk immer- 
sion. In Curir (2006) we showed that our embedding 
procedure is not triggering, per-se, a bar instability, and 
as well, no instability is triggered because of our chosen 
numerical resolution. 

The units of our simulations are: lO^'^h"^ Mq for 
the mass, Ih^^kpc for the length, 0.98h^^Gyrs for the 
time. 

2.1 The baryonic disk 

The spatial distribution of particles follows the expo- 
nential surface density law: p — poexp— (r/ro) where 
iQ is the disk scale length, rg = Ah~^ Kpc, and po 
the surface central density. The disk is truncated at 
five scale lengths with a radius: Kdisk =20h~^ Kpc. Its 
mass is 1.9 x 10^°h-'^MQ. 

This value of mass has been chosen and employed in 
the previous papers as the critical one for the bar in- 
stability: the disk/halo mass ratio for this case is in 
the instability range from the classical point of view, 
the presence of a cold gas knot originated by gas cool- 
ing does stabilize it against the bar formation, and it 
becomes unstable again if the gas is allowed to form 
stars. 

To obtain each disk particle's position according to 
the assumed density dis t ributi on, we used the rejec- 



tion method iPress et al.l (|1986[ ) . We used 56000 star 
particles and 56000 gaseous particles to describe our 
disk. The masses of gas, stellar and DM particles are 
respectively 7 x lO'' /i-^M©, 2.8 x IO^/i-^Mq, 1.194 x 
10^h~^ Mq. The (Plummer equivalent) softening 
length, the same for DM, gas, and star particles and 
black hole, is 0.5/i~^ kpc. 

In the cosmological case we embed the disk in the 
high-resolution cosmological simulation, at redshift 2, 
in a plane perpendicular to the angular momentum vec- 
tor of the halo and in gravitational equilibrium with the 
potential. Its center of mass corresponds to the posi- 
tion of the DM particle having the minimum value of 
gravitational potential. The initial redshift corresponds 
to 10.24 Gyrs down to z = in our chosen cosmology. 
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During the evolution, new star particles are formed 
from the gaseous particles. We refer to such new 
component as 'new star component', whereas we will 
call 'old star component' the non dissipative particles 
present in the disks at z = 2 (at t = in the isolated 
cases). 

In the disk center of mass resides a BH which is rep- 
resented as a coUisionless, heavy mass point. It is able 
to swallow the gas in its neighborhood. The feedback 
parametere/ assigns a certain percentage of the accre- 
tion energy to be thermally coupled with the gaseous 
component. This feedback regulates the BH accretion, 
which, without such effect, is catastrophic. The default 
value of such parameter is ef — 0.05 (see below), but 
we will explore in the paper other diff erent values. We 



used the BH/energy feedback model of lDi Matteo et al 
This model hypothesizes that the gas accretion 
is limited by the Eddington rate: 



M = 



AnGMBHrrip 



where Mbh is the BH mass, the proton mass, 
(Tt the Thompson cr oss section and t r the radiative 
efficiency. We follow iDi Matteo et al. I (|2003) and set 
e,. = 0.1. A fraction e/ of the energy emitted by the 
BH couples with the surrounding gas, giving an energy 
output rate of Ef^ed = ^f^rMsHC^- The parameter e/ 
is the feedback parameter. 

2.2 Star formation recipe 

We use the sub-grid star formation prescription by 
Springel & Hernquist (2003). In it, when a gas particle 
overcomes a given density threshold its gas content, is 
considered to reside in a multi-phase state, that cor- 
responds to the equilibrium solution of an analytical 
model describing the physics of the multi-phase inter- 
stellar medium. Such a solution gives an effective tem- 
perature for the gas particle, obtained as a weighted av- 
erage of the (fixed) temperature of the cold phase and 
of the hot phase. The temperature of the hot phase 
is set by supernova feedback and by the efficiency of 
cloud evaporation. The resulting inter stellar medium 
(ISM) has an effective equation of state Pe//(p), which 
is stiffer than isothermal and prevents Toomre instabil- 
ities even when a large amount of gas is present. When 
this prescription is used, the star formation rate agrees 
with the Schmidt law in disk galaxies, as obtained e.g. 
by Kennicut (1998). The effective temperature drives 
hydrodynamical interactions of the gas particles. More- 
over, the model consistently gives a star formation rate, 
which is used to spawn a star particle from the gas 
one on the basis of a stochastic prescription; the initial 



mass function adopted is the Salpeter's one. The star- 
forming from the gas has the same position and velocity 
as the gaseous particle. 



3 Simulations 

We performed 11 simulations of a disk-|-halo isolated 
system, and 2 Cosmological simulations, where the disk 
was embedded in our selected cosmological halo. We ex- 
ploited the non-public versi on of the parall el Tree+SPH 
N-body code GADGET-2 (ISpringellboOsI) (courtesy of 
V. Springel). The simulations run on the CLX com- 
puters located at the CINECA computing center (BO, 
Italy ). The configurations have been evolved for 10.44 
Gyrs (tfin). 

The main parameters and the final properties of our 
set of simulations are listed in Table [1] We kept fixed 
the disk mass and the gas to stellar mass ratio, and 
varied the initial BH mass and the feedback strength 
€f. We focused on the effect of the presence of a 
BH on the bar instability, and re f er the reader to ou r 
previous wo rks Curir et al. ( 20061 ). Curir et al. ( 2007 ). 
Curir et al.l (j2008l ) for a study of the impact of other 
significant parameters of our disk-|-halo system on such 
instability. 

As a geometrical global measure of the bar strength, 
we defined the value of the ellipticity, e — 1 ~ b/a, 
(where a, b are the semi-axes of the contour) chosen on 
the iso-density level where such a value is maximum. 

We defined the bulge as the total mass contained 
inside a radius centered on the new stars center of mass 
and corresponding to a minimum thickness of Ah~^kpc. 
The bulge s in our simulations ar e "pseudobulges" as 
defined in iKormendv et al.l ( 2002 ) . The initial disk is 



bulgeless and the bulge is formed by secular instability 
and accreted by the formation of the new stars. In 
our simulations the timescale of the formation of these 
pseudobulges is 3-4 Gyrs. The final masses consisting 
of old and new stars are between 6 x 10^h~^MQ and 
7 X IQ^H-^Mq. 

In Table [1] we present the simulation number (I col- 
umn), the initial mass of the BH in code units (II col.), 
the feedback value, i.e. the percentage of accretion en- 
ergy coupled with the gas heating (HI col.). Moreover, 
as final values we present the maximum ellipticity for 
the old star system(IV col.) and for the new star (V 
col.), the final mass of the accreted BH (VI col.), the 
major axis corresponding to the maximum ellipticity for 
the old stars (VII col.)and for the new stars (VIII col.), 
the central mass concentration CMC defined as the 
mass of all the components (stars, new stars and gas) 
inside a radius of 1.5h~^ kpc followed by the mass frac- 
tion of BH to the CMC (IX col.) and the time needed 
for the accreting BH to reach the value lQ^h~^MQ . 
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Table 1 Simulations: o.s.=oId stars , n.s.= new stars. Masses are in units of 10^''^ ^Mq, lengths in h ^ kpc, times in 



Gyrs 


N 


BH rriini 






e (o.s.) 


e (n.s.) 


BH m fin 


^max (O-S.) 




CMC — BHfrac 




11 


10-6 







0.747 


0.75 


0.18 


6.37 


1.6 


0.23 - 0.72 


6.289 


12 


0.001 







0.43 


0.44 


0.23 


6.28 


6.27 


0.3 - 0.76 


2.05 


13 


0.005 




10-4 


0.845 


0.85 


0.021 


7.7 


1.2 


0.19 - 0.11 


0.993 


14 


lo-'^ 




0.05 


0.8 


0.81 


6.8 X 10-5 


7.94 


1.5 


0.22 - 10-4 


i ^ tfin 


15 


10-4 




0.05 


0.69 


0.71 


1.15 X 10-4 


6.16 


6.6 


0.14 - 8.2 X 10-4 


i ^ if 171 


16 


10-4 




0.01 


0.74 


0.82 


2. X 10-4 


3.45 


1.31 


0.15 - 1.3 X 10-'"^ 


i ^ if 171 


17 


0.001 




10-5 


0.58 


0.49 


0.21 


6.59 


3.73 


0.16 - 1.3 


0.1667 


18 


0.001 


5 


X 10-5 


0.77 


0.69 


0.027 


8.29 


2.49 


0.19 - 1.42 


1.389 


19 


0.001 




10-4 


0.73 


0.75 


0.014 


8.16 


1.96 


0.18 - 0.078 


5.49 


110 


0.001 




0.001 


0.83 


0.83 


0.0024 


6.72 


1.21 


0.16 - 0.015 


i ^ if 171 


111 


0.001 




0.01 


0.74 


0.85 


0.0018 


7.8 


1.36 


0.17 - 0.01 


t > if in 


cl 


0.001 


5 


X 10-5 


0.53 


0.477 


0.0065 


8.18 


4.08 


0.15 - 0.04 


i > if in 


c2 










0.476 


0.48 





7.2 


3.8 


0.1 


i > if in 



4 Results 

Our main result Is that In all of our simulations, a stel- 
lar bar Is still living at the final time In the old star 
component. The new star component at the end of our 
runs Is arranged In a bulge component which presents a 
more or less barred shape depending on the Initial mass 
of the BH and on the feedback. Overall, the presence 
of a BH Is never able to quench the stellar bar, not 
even when Its mass Is significantly larger than allowed 
by observations. 

An example of the final configuration Is shown In 
Fig. [T] for simulation 18. Fig [T] has been constructed 
with the same box-size, number of levels, and density 
contrast as in Papers 1, 2 and 3. 

In Fig. 2 we show the value of the elllptlclty of the 
bar at the final time, as function of the final mass of 
the BH. Only results for simulations having the same 
initial masses of the BH {10^ Mq) are shown here. 
The final BH mass Is a direct function of the feedback 
parameter, the larger one corresponding to a simulation 
with no BH feedback. This diagram show a clear, if ob- 
vious, correlation for the masses larger than IO^/j-^Mq: 
more massive black hole-weaker bar. In particular, the 
point of lower elllptlclty in such figure correspond to 
the final time of simulation 12, where the feedback pa- 
rameter is zero and the BH mass is 2.3 x 10^h~^ Mq, 
a value higher than the ni lnlmum mass requested by 
Hozuml fc Hernguistl (|2005l) for a complete bar destruc- 
tion. In our case the bar instability is not totally inhib- 
ited, but the elllptlclty is the weakest one. 

In Fig. 2 the peak in the elllptlclty corresponding 
to a feedback ef — 0.001 (simulation 110) is created by 



the competition between two effects: the star forma- 
tion which is working in favor of the bar (and which 
is weakened by the BH feedback) and the increase of 
the central BH mass -I- CMC, which is inhibiting the 
bar. The star formation rate Increases moving from the 
first to the second point of diagram 2 (because the BH 
feedback decreases), it is stationary going to the third 
point (but the central BH mass -I- CMC has Increased 
enough to weaken the elllptlclty) and it decreases again 
after. 

More details about the Impact of the BH on the el- 
llptlclty are in Fig. 3, where we see the final elllptlclty 
profile, as a function of the semi-major axes of the Iso- 
denslty contours for the simulations of Fig. 2. 

In Table 1, we see that the values of the CMC are 
different for different simulations and in most of the 
cases they are bigger, or much bigger than the BH fi- 
nal mass. Therefore one can deduce that the effect of 
the black hole in the disk center is not really a direct 
gravitational Impact on the bar, but more an effect of 
modulating and, for the more massive BH and lower 
feedback values, enhancing the CMC formation, which, 
in turn, has an effect on the bar elllptlclty. In Fig. 4 
we have shown the direct correlation between the final 
BH mass and the CMC: the more the black hole is mas- 
sive, the higher is the CMC. In Fig. 4 we present only 
the cases starting with the same initial BH mass: as in 
Fig. 2, more massive BHs correspond to lower values 
of the feedback. In the same figure, we show an hori- 
zontal line: the value of the CMC for the case of the 
disk without BH. It is evident that for the higher values 
of the feedback (first three points of the diagram), the 
value of the CMC is reduced by the gas heating, when 
compared with the CMC of the same disk without a 
central BH. The difference in CMC between cases with 
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Fig. 1 iso-density contours of old stars in the planes xy, xz and yz (left three panels) and new stars (right panels) for 
simulation 18, at the end of the evolution 




Fig. 2 Ellipticity as function of the final mass (in solar masses) of the BH for the simulations with the same BH's initial 
mass (W'^h-^MQ). 



strong feedback and weak one is almost completely due 
o the new stars component and to the gas content. 

We also compared the effect of having different evo- 
lutionary patterns of our central BH on the onset of 
the bar, i.e. different initial BH masses with different 
feedback efficiencies. With this aim, we looked for the 
time t at which the mass of the BH is equal to the value 
lO^/i"^M0, independently of the initial BH mass. We 
then plotted in Fig 5 the ellipticity as a function of the 
time t needed for the BH to reach such a mass value. 

We can disentangle the effect of the bar evolution by 
comparing the bar ellipticities measured in the same 
disk, at the same times, in a simulation without a cen- 
tral BH (dotted line in Fig. 5). 

This comparison shows that a too fast or too slow 
growth of the BH results in a small effect on the bar 
ellipticity. A significant change of ellipticity, with re- 
spect to the simulation without the BH, only happens 
if the BH takes from 2 to 5 Gyrs to reach a mass of 
10® Mq. Of course, the subsequent evolution of the bar 
does depend on the mass of the BH; a BH which quickly 
reaches our reference mass and then continues to grow, 
will have a larger effect on the bar ellipticity at the final 
time, as shown in Figure 2. 

In Fig. 6 and 7 we provide the profiles of the differ- 
ences of the ellipticities with and without the BH, for 
the old (Fig. 6) and the new (Fig. 7) stars at the same 



times as in Fig. 5., i.e. when the BH mass reaches the 
value 10® Mq. These differences are always negatives, 
for the new stars except for t=2.05 (case 12 in Tab. 1), 
meaning that for all the other times, the bar is stronger 
at all radii if the BH is not present. 

For the old stars profiles the trend is less evident: at 
some times we have radii in the disk where the elliptici- 
ties of the bar containing the BH is slightly higher. But 
it is evident that the differences are always negative at 
all times as far as the internal radii is concerned. We 
interpret this as the confirm that the effect of the BH is 
to decrease the strength of t he inner bars. This is con 



sistent with the findings of iHollev-Bockelmann et al 



( 2002f ) which however only considered the stellar com- 
ponent in a triaxial system and used an adiabatically 
growing potential to simulate the presence and growth 
of the central BH. They noticed that the major impact 
of the BH on the triaxiality of the galaxy involves at 
maximum the 10% inner mass of the galaxy. We esti- 
mated indeed that the major impact on the variation of 
the ellipticity is on the inner iso-density contour, con- 
taining Rs 10% of the mass, when the BH is present 
with respect to our control case. 

These results are qualitatively confirmed even if we 
evaluate the ellipticity using the inertial tensor inside 
our stellar system. 
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t(Gyr) 



t(GyT) 



Fig. 5 EUipticity of the bar when the accreting BH reaches the value l(f Mq , as function of the time accretion (solid 
line). EUipticity at the same time for a disk without central BH (dashed line). The ellipticity is evaluated for old star (left 
panel) and new stars (right panel) 



CMC versus BH 



ellipticity 




Jb-1.E-2 fb-f.e-3 fb-1.e-4 fb-5.E-5 fb-1.e-5 fb-0 ^ 



B to 12 



ellipticity 




fb^l B-2 Jb,l.s-S Jb, 



fb=3 e-5 fb=1.E-S fb = 



Fig. 3 Profiles of the ellipticity at the final time for the 
set of simulations with the same initial BH mass 10^ h"^ Mq 
and different feedback. The profiles are for the old stars 
component (upper panel) and for the new one (lower panel) 



Fig. 4 Final central mass concentration (in solar masses) 
versus final BH mass (in solar masses), for the cases having 
the same initial BH mass. The horizontal line indicates the 
value of the final CMC in the case of the disk without BH. 




rmmajor axis 



Fig. 6 Differences between ellipticity profiles of the stellar 
bar (old stars) with (e(BH)) and without (e) the BH, when 
the accreting BH reach the value lO^/i'^M©, the different 
colors correspond to different accretion times. 
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4.1 Cosmological cases 
4- 1.1 bar strength 

We performed two cosmological simulations, one with 
a BH having a mass Mbh = lO'^ Mq (cl) and an- 
other without the BH (c2). We chose ej — 5.10~^ as 



our feedback param eter 



cedure described in Curiretal 



We us e d the embedding pro- 
(|2006h . The values of 



the ellipticities show a very mild impact of the BH on 
the onset of the bar instability and its evolution. In 
Table 1 we reported the final (z — 0) value of the bar 
ellipticity, which is very similar in the two cosmological 
cases both in the old and in the new stellar component. 
Note that, for the cosmological run, we used a quite 
high BH initial mass and a low feedback parameter. 
Such parameter choice is the one which shows to be the 
most effective in dumping the bar instability in our iso- 
lated cases (case 18 in Table 1). Even so, the presence 
of the central BH has practically no inhibiting effect on 
the bar instability. 

4-.1.8 DM concentration 

In this section we investigate if the BH has an impact 
on the concentration parameter of the haloes. We then 
compare three different haloes: the halo from our orig- 
inal DM-only simulation, selected for the disk immer- 
sion; the same halo containing the disk, and the same 
halo containing the disk and the BH. 

In the table [2] we show the behavior of the NFW 
concentration parameter as a function of the redshift 
z in the three cases. It is known that the baryonic 
matter has an influence on the halo concentrati on in- 
creasing it by an amount that iLin et al.l (|2006[ ). e.g, 
estimated to be around 10%. Moreover, the concentra- 
tion evolves approximately linearly with redshift (see 

new stars 




e.g. Bul lock et al.l (j200ir) l: at redshift z = 2 we expect 
a concentration parameter of the order of one third that 
found at z = 0. Our DM halo shows a concentration 
evolution which is evolving more than linearly: it in- 
creases by a factor « 4 from z = 2 to z = 0. In the 
Table, we also show the increase in the concentration 
of the DM halo due to the presence of the baryonic 
disk only (case cl). This increase is higher than the 
expected 10%. On the other hand, from the table it is 
clear that the presence of a BH has practically no effect 
on the concentration parameter. 



5 Kormendy relation 

Figure [8] shows the values of the stellar velocity disper- 
sion inside the bulges at the end of the evolution as 
function of the BH masses. 

We co mpared data from our sim ulations with obser- 
vations bv lKormendv et al.l (j2002r) . Data from our sim- 
ulations are located in the pseudo-bulges region, with is 
consistent with the origin of the bulge structure in our 
simulations: they are not classical hot bulges, rather, 
they derive from secular evolution of the bar. 

We can't find the phenomenologically observed cor- 
relation of velocity dispersions with the BH masses. It 
must be noticed that we vary the BH mass and feed- 
back parameter within the same DM halo, so from this 
point of view our Fig. 8 simply suggest that varying 
the feedback efficiency and the initial BH mass has no 
effect on the stellar velocity dispersion, which remains 
a reasonable proxy for the halo mass. 

The values that better agree with observation are 
obtained from simulations using a high BH feedback 
parameter. This trend suggests an useful constraint for 
such parameter in the simulations. Another possibil- 
ity is that pseudo-bulges could behave differently, in 
this regard, from hot bulges. This point remains to be 
deepened both observationally and theoretically. 



Fig. 7 Same as in the previous figure, but for the new 
stars. 



6 Discussion and conclusions 

We have presented 1 1 isolated simulations and 2 cosmo- 
logical simulations of the evolution of a stellar -|-gaseous 
disk embedded in a DM halo,with the same disk-to- 
halo mass ratios. The aim was to evaluate the impact 
of a massive black hole, with different accretion histo- 
ries on the formation and evolution of a stellar bar. 

The old star component shows a long-lasting bar, 
10 Gyrs old, in all our simulations, regardless of the 
presence, the mass and the energy feedback efficiency 
of a central BH. 
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Table 2 Halo concentrations 



z 


halo 


cl 


c2 


2 


4.5 


4.5 


4.5 


1 


8.8 


12.8 


12.9 





19.7 


29.1 


29.9 



We noticed a mild impact of the BHs having masses 
greater than IO^H^^Mq on the old star component cl- 
lipticity, impact which appear stronger on t he bar ellip- 
ticity of the new stars. In a previous paper (jCurir et al 
( 20081 )) we found that the star formation, by reduc- 
ing the central gaseous mass concentration, allows the 
bar to survive until the end of the e volution in m assive 
disks, at variance with the results inlCurir et al 



mm 

( where the gas was not allowed to form stars): in the 
latter case a gas fraction 0.2 was able to destroy the 
bar. 

In this work we show that a BH placed and growing 
at the center of the disk has a small influence on the 
ellipticity of the bar developing in the preexistent star 
population. On the other hand, the BH with its feed- 
back has an action on the formation of the CMC, which 
has a major role in quenching the bar. The bar formed 
by the newly formed stars is always weaker than in the 
disk without BH, independently on the BH's evolution 
and feedback. Also in this case, however, the presence 
of the BH does not completely stop the bar formation. 

A direct comparison between t he results of our TV- 
body s imulations and the results of Hozumi fc Hernguist^ 
( 20051 ) IS not possible. The reason is that their numer- 
ical model consists in a razor-thin disk model, without 
bulge an d halo. It is evolved us i ng a s elf-consistent field 
method ( Hernguist fc Ostriker ( 19921 )). and not with a 
A'^-body code. Whereas the Hozumi's disk is stabilized 
by a black hole having mass of IO^'^Mq our disk forms 
a bar even if the central BH has a higher mass. 

Adding gas and star formation to the same disk, we 
notice that an accreting massive BH has an impact on 
the bar ellipticity but it is not able to destroy it com- 
pletely. 

The effect of the BH is more evident on the inner 
bar ellipticity, which is reduced also when the effect of 
a star-forming gaseous component is considered. 

A comparison with the phenomenology confirms that 
the default feedback parameter value 0.05 seems the 
more appropriate to reproduce the correlations between 
the velocity dispersion and the black holes masses ob- 
served by Kormendy. 
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1.8 2.0 2.2 2.4 2.6 

log CTg {km s"'') 

Fig. 8 The BH masses versus the velocity dispersions 
in our simulations (green points) compared with the Ko- 
rmendy observed data (cit) 
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